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We have studied the dynamics of a cylindrical column of anisotropic, charged fluid which is 
experiencing dissipation in the form of heat flow, free-streaming radiation, and shearing viscosity, 
undergoing gravitational collapse. We calculate the Einstein-Maxwell field equations and, using the 
Darmois junction conditions, match the interior non-static cylindrically symmetric space-time with 
the exterior anisotropic, charged, cylindrically symmetric space-time. The behavior of the density, 
pressure and luminosity of the collapsing matter has been analyzed. From the dynamical equations, 
the effect of charge and dissipative quantities over the cylindrical collapse are studied. Finally, we 
have presented the solutions for this dissipative, cylindrical gravitational collapse and have discussed 
the significance from a physical standpoint. 



I. INTRODUCTION 



An important task in relativistic gravity and astrophysics is to build realistic models to describe the final fate of 
a star which has exhausted its nuclear fuel GUI]- While the star was still burning, the material was held up by the 
radiation pressure generated inside its core [3|. Once the fuel is exhausted, the star begins its collapse under the 
internal pull of gravity. The fate of any gravitational collapse depends upon the mass of the collapsing star A 
star with a mass of the order of the Sun will eventually collapse into a white dwarf, where the material is held up 
bTj' by the electron degeneracy pressure generated by electrons moving within it at high speed. A star which has about 
l— "™ ', four to five times the solar mass will however settle into the state of a neutron star, in which even the atoms have 
■ collapsed under the enormous pressure of gravity, leading to the formation of minute electrically neutral regions inside 
it, of dimensions greater than the size of an atom [5J. This object made up purely of neutrons will be held up due 
to the neutron degeneracy pressure and will be about twenty kilometers wide. Even more massive stars will neither 
settle into a white dwarf nor a neutron star and the gravitational pull will be sufficient to collapse it into a black hole. 
These are one of the most exotic objects in the universe as they collapses to within the " Schwarzschild radius" and 
' therefore no information is able to escape from within, what is called, the event horizon. 

There has been an extensive study of collapse of dust and fluid under gravity since the formulation of General 
Relativity by Einstein. Contributions by S. Chandrasekhar, who worked out the maximum mass of a stable white 
dwarf star, known as the Chandrasekhar limit, to the collapse study of adiabatically flowing dust by Oppenheimer 
and Snyder Q , set the standard for research in this area. P. C. Vaidya studied the external gravitational field of a 
£> | stellar body giving out radiations. Misner and Sharp M, [101] studied spherically symmetric collapse and contributed 
by giving an expression for the total energy. Others [ill ] studied the geometry of the Vaidya metric. K. Lake, L. 
Herrera, N.O. Santos and many others [12h21| have studied different cases of spherically symmetric fluids undergoing 
collapse. A realistic model of such a collapse should include a radial heat flux as the process of gravitational collapse 
is highly dissipative [22], [22j ■ Recently, Sharma and Tikekar [24[ has studied the evolution of non-adiabatic collapse of 
a shear-free spherically symmetric stellar configuration with anisotropic stresses accompanied with radial heat flux. 
The effects of electromagnetic field on the energy density inhomogeneity in the relativistic self-gravitating fluids for 



spherically symmetric spacetime has been studied by Sharif and Bashir 25| . There has also been extensive studies on 
quasi-spherical and self-similar collapse [26| over the years. An extensive list of papers on the study of gravitational 
collapse can be found in [l[ and (24|. 

The study of non-spherical gravitational collapse has gain ed in momentum following the discovery of cylindrical 
and plane gravitational waves. Bronnikov and Kovalchuk [27] were the pioneers in the investigations on non-spherical 
gravitational collapse. Ashtekar and Varadarajan obtained the expression for the energy per-unit length (along the 
symmetry direction) of gravitational waves with a spacelike symmetry in 3+1 dimensions, with cylindrical waves as 
a special case [28[ . The definition of cylindrical symmetry in general relativity as a specialization of the case of axial 



symmetry has been considered by Carot et al. [29( and Barnes [30(. Chiba [31| studied the case of cylindrical dust 
collapse. Others |32H34j investigated various aspects of cylindrical collapse of counter rotating dust and rotating 
cylindrical shells. Hayward studied gravitational waves, black holes and cosmic strings in cylindrical symmetry. 



Considering the most general vacuum cylindrical spacetimes, Goncalves 35[ has demonstrated that there are no 



trapped cylinders in the spacetime, and has presented a formal derivation of Thome's C-energy, based on a Hamiltonian 
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reduction approach. Di Prisco et al. |3jJ studied shear-free cylindrical gravitational collapse for an interior non- 
rotating fluid with anisotropic pressures and exterior vacuum Einstein-Rosen spacetime. The case of the collapse of a 
heat conducting charged anisotropic fluid cylinder have been studied by Sharif and Abbas [38| and that of a charged 
fluid cylinder with shear viscosity by Sharif and Fatima [39| . 

In this paper, we have extended the work done by Sharif and Abbas [38j j and Sharif and Fatima (39j and have 
examined the effect of charge, heat flow, radiation and shear viscosity on the gravitational collapse of a cylindrical 
column of anisotropic fluid. We have formulated the field equations for the source and have employed the Darmois 
junction condition for the smooth matching of the interior non-static cylindrically symmetric space-time and the 
exterior anisotropic, viscous, charged, cylindrically symmetric space-time at the boundary S, to derive a few results. 
The paper is organized as follows: In the next section, we describe the gravitational source along with the corresponding 
physical parameters like the expansion scalar, acceleration, shear tensor and the Einstein-Maxwell field equations. In 
section III, the exterior metric and the junction conditions are given. The results, including the dynamical equations 
and the solutions are presented in section IV. The summary of this whole exercise is given in section V. 

II. THE INTERIOR METRIC AND THE FIELD EQUATIONS 

We consider a collapsing cylinder filled with an anisotropic, charged fluid and undergoing dissipation in the form of 
heat flow, free-streaming radiation, and shearing viscosity, bounded by a timelike cylindrical three-surface S, which 
divides the space-time into two distinct 4-dimensional manifolds V + and V~ . 

A. The interior spacetime 

For the interior V~ space-time we take the general non-static cylindrically symmetric metric in the comoving 
coordinates given by 

ds 2 _ = -A 2 (t,r)dt 2 + B 2 (t,r)dr 2 +C 2 (t,r)(d9 2 +dz 2 ) (1) 
In order to represent cylindrical symmetry, the range of coordinates is required to be as follows: 

— oo < t < oo, < r < oo, 
< 9 < 2tt, -oo < z < +oo, 

with the coordinate labels, x° = t, x 1 — r, x 2 = 9 and x 3 = z. The interior energy-momentum tensor is given, 
according to relativistic hydrodynamics, as 

T a p = (n + P±)V a Vp + P±g a p + (P r - P±)XaXp + VaQp + Vpq a + eljp - 2i]a al3 , (2) 

where, fi — ¥ energy density, P± — > tangential pressure, P r — > radial pressure, q a — > heat flux, e — > the radiation 
density, V a — > 4-velocity of the fluid, x° un it 4-velocity in the radial direction, l a —> a null 4-vector and r/ — > 
coefficient of shearing viscosity > respectively. 

The shear tensor a a p, the 4- acceleration a a and the expansion G are defined as 

<y a! 3 = ^ ((K^ + Vp ;a ) + (flaVa + dpV a )) ~ ^{g a fj + V a Vp), (3) 



a a = V aiP V p , (4) 



and 



6 = V£. (5) 
Since we have assumed comoving coordinates for the interior metric, we have 

V a = A- 1 S%, x a = B- 1 5«, l a = A- x 6% + B- 1 5%, q a =B~ l q5^, (6) 

such that 

v a v a = -i, x a x a = h x a v a = o, q a V a = 0, l a V a = -l, l a l a = Q. (7) 
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In view of the equations Q, ([5]) and ©, we obtain the acceleration and the expansion scalar as follows: 

(8) 

Using ^ to ©, we obtain the non-zero components of the shear tensor as 

2 1 

(Til = -J=B 2 (T 7 (T 2 2 = 033 = -— C 2 a. (10) 




The shear scalar a defined by 



is therefore obtained as 



a 2 = layer* (11) 




(12) 

As defined by Chiba [3l| . a cylindrically symmetric spacetime may be defined locally by the existence of two 
commuting, spacelike, Killing vectors, such that the orthogonal space is integrable. For such a spacetime, there exist 
coordinates (9,z) such that the Killing vectors are (£e,£ z ) == (d/d6,d/dz). The existence of cylindrical symmetry 
about an axis implies that the orbits of one of these vectors are closed but those of the other is open. Each of these 
Killing vectors must be hypersurface orthogonal. The norms of these Killing vectors are invariants |4p| . namely the 
circumferential radius 



and the specific length 



with £( 2 ) = 9g, £(3) = d z , under the sign convention that spatial metrics are positive definite, the dot representing 
contraction and the flat b represents the covariant dual with respect to the space-time metric. The gravitational 
energy per specific length in a cylindrically symmetric system (also known as C-energy) as defined by Thorne |4lj 
and modified by him to render it finite in spacetime, is given by 

E= ~(l-r 2 V<W f), (13) 

o 

for which 

r = (l, 

where the areal radius is f and E is the gravitational energy per unit specific length of the cylinder. 

Analogous to the Misner and Sharp energy for spherical symmetry [42j . the specific energy of the cylinder due to 
the electromagnetic field is therefore given by 




B. Electromagnetic energy tensor and Maxwell's equations 

The electromagnetic energy-momentum tensor for the charged fluid is given by 

rS" 5 = h fa* - \ F1&F ^() ■ (is) 
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and the corresponding Maxwell's equations are 

F af 3 = 1pp, a ~ 1p a ,p, (16) 

F al5 . = 4nJ a , (17) 

where F a p is the electromagnetic field tensor, ip a is the corresponding four potential and J a is the four current 
density vector. Since the charge is comoving with the fluid, there is no magnetic field, and the four current density is 
proportional to the four velocity i.e. we have 

V>a=#S, J a = P V a , (18) 

where rp(t,r) is an arbitrary function and p(t,r) is the charge density. 
From the Maxwell's equations we obtain 

A' B' _C'\ ,, . .„ 2 



V/' - I -j + — - 2-1 t/y = AirpAB , (19) 

V/-(| + f-2§)v/ = 0, (20) 

where the first equation is for a = and the second is for a = 1. Here the dot and the prime represent the partial 
derivatives with respect to t and r respectively. Integrating (|19[) we obtain 

,, 2s AB , . 

V = -03-. ( 21 ) 

where 

s(r) = 2tt / pBC 2 dr (22) 
Jo 



is the total charge distributed per unit length of the cylinder. Equation (|21[) is in conformity with the conservation 
of charge and satisfies Eq. (|2"0")) . 

C. The Field Equations 

We now find the field equations for this distribution of fluid. The Einstein field equations for the interior metric 
can be written as 

G-^MT-p + T^-) (23) 

where G~p is the Einstein tensor for the interior metric. There are five non-zero components of (|23[) for the metric 
([1]) with energy- momentum tensor $Z§, which are 

Goo = 8^(r oo +T ( emr ) 



i.e. 



4s 2 A 2 C ( B C\ ( A\ 2 ( C" C ( B' C 



8 ^ + ^ A + — =C[ 2 B + C + {b {- 2 - + c( 2 B-c"- (24) 



Similarly, 



which yields 



(em) 



M9 + C)AB = 2 . (25) 
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The remaining equations are 

G u = 87r(T-+T/r r ) 



4 \ „, 4s 2 B 2 , . , 4s 2 B 2 



8vr \P r + e - -= V a B 2 - — - = Sn(P rcff + e)B 



B\ 2 (C ( cV iC*\ {C'V A'C 



= -<I 2 c + c - 2 lc + c +2 ac' (26) 



where the effective radial pressure is defined as 



4 

Pr eft = Pr-^W 



and 



G 2 - 2 = 8^(T 2 - 2 + Tj 2 eror ) = Stt (Pj. + c 2 + |J = 8nP^ ff C 2 + |J 



C*\ 2 / B C A IB C\ BC 



A) \B C A \ B C BC 
C\ 2 fA" C" 



B ) \A C A 
with the effective tangential pressure as 





a 


\ B'C'\ 


(1 


c 


J BC ) 




p± + 


2 



(27) 



III. EXTERIOR METRIC AND THE JUNCTION CONDITIONS 

Exterior to the hypersurface E in the 4D manifold V + , we consider Vaidya's metric Q in presence of charge in 
the retarded time coordinate as considered by Chao-Guang (43|, but with a signature flip. The introduction of the 
retarded time coordinate removes the singularities of the original line element. Let M(u) and Q(u) be the mass and 
charge of the fluid respectively inside the hypersurface E, where u is the retarded time coordinate. Then the exterior 
field in this cylindrically symmetric spacetime can be defined as 

da 2 . = -( ~ 2M(M) + 9^1) du i _ 2dRdu + R 2 {d9 2 + dz 2 ). (28) 

The intrinsic metric for the hypersurface E which enables a description in comoving coordinates of the interior 
spacetime, is given by 10] 

(ds 2 ) s = -dr 2 + Rl(r){d9 2 +dz 2 ), (29) 
where ()s means the value of () on E and £ z = (r, 9, z) represents the coordinates on E, i.e. 

To match the interior and the exterior space-time, we follow the prescription of Darmois and Israel pij which 
demands: 

• The first fundamental form must be continuous over the hypersurface E i.e., the continuity of the metrics as 
V ± approaches E: 

(cfe 2 ) s = (ds 2 _h = (d«+)E. (30) 

• The continuity of the second fundamental form. This gives the continuity of the extrinsic curvature Kij at the 
hypersurface E: 

= K+ - Kr = 0. (31) 
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According to Eisenhart |45| , the extrinsic curvature of E is given by 



Kt = -n 



(<7,H,v = 0,1,2,3). 



(32) 



where are the outward unit normal vectors to the hypersurface E, x ±A1 are the coordinates of V^. 

On the r = constant hypersurface, we have dr = 0. Using this condition in ([1]) and comparing with (|29j) keeping in 
mind the junction condition (|3T))) . we get, 



d7 



i? s (r) = C(t,r E ). 



(33) 



We may also write the exterior metric (|2"8"j) as, 



(d4)s = - 



-2M{u) Q 2 {u) 



2dRj 
du 



du 1 



Rl(d8 2 +dz 2 



Now, using the junction condition (|30l) and matching with the metric on the hypersurface E, we get 



du 
dr 



-2M(u) 
Ry, 



2 (u) 2dR 1 



-,-1/2 



Rl 



du 



(34) 



(35) 



In order that we can apply the junction conditions, we require that E has the same parametrisation whether it is 
considered as embedded in V + or in V~. In the coordinates of the interior spacetime V~ , the bounding surface E 
will have the equation 



f(t,r) = r - r s = 0, 

where is a constant. 

Since the vector df/dxl is orthogonal to E, so the unit normal vector to E in the X- coordinate system is, 

n- =0,B(t,r E ),0,0. 
In the coordinate system of V + , the equation for the surface E may be written as, 

f(u,R) = R-Rj;(u) = Q. 
The vector df / dx% , orthogonal to the hypersurface E is therefore given by, 



Of 
dx+ 



dR^ 
du 



1,0,0 



So the unit normal to E in the V + coordinate system is, 



-2M{u) 



l {u) 2dR Y 



du 



-1/2 



dRy 
du 



,1,0,0 



The extrinsic curvature for the hpersurface E in the V + coordinates as calculated by using (|32[) is given by 

KZ 



-( — 
\AB 



K, 



22 



K. 



33 



S 



K, 



no 



K+ = 



d 2 u ( du" 1 
aW 1 KaW, 



R2~R3J \~dr~) 



R™ + (91- 2M ) ^ 
dr \ R ) dr 



(36) 

(37) 
(38) 

(39) 

(40) 

(41) 
(42) 

(43) 

(44) 
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On account of the continuity of the second fundamental form given by (|31[) . we obtain the following relations on 
matching @TJ to (O and (O to (gU) 



cPu 
ah 2 



M 
R? 



R 



dR 



9l 

91 

R 



- 2M 



du 
ah 



AB 
f CC 



(45) 
(46) 



IV. RESULTS 



dr J \ i?s J \dr J \ dr 

Putting this value in (l4"6l and using ([55]). we have 



We now use the relations obtained above and simplify them to find useful results. From (|35j) we have by rearranging, 

^)(^-2m)=rJ^Y 1 -2rJ^). (47) 



(48) 



du\~ 1 c a 



dr J \A B I 

Again using (J33]) and squaring (|48[) we obtain the total energy entrapped inside the surface E as follows: 




2 „ x 

C"\ 2 \ Q 2 



BJ J ' 2C- (49) 

Taking the interior and exterior charge to be the same on the hypersurface S (i.e. Q = s) and using (j49|) and (fM)) . 
we obtain 

E'= l -+M, (50) 

which indicates that the difference between the two masses is equal to 1/8, as obtained in [38[ and (39|, which is a 
consequence of the least unsatisfactory definition of C-energy due to Thorne [4l[ . 
Using the expressions (|8|), ((9]) and (fT2|) we can reconstruct (|25|) as follows: 

4^(g + e) = i Q(9 - V3a)' - VSa^-) . (51) 



Differentiating (|48|) with respect to t and substituting in (j45|) . we obtain the following expression with the help of 
([38]l and (01), 

' ^ Ui ,52, 





i- 






f C" 


BC'\ 


1 


a 2 








I c 


BC ) 


~ AB 



A B C 2 \C 

Using (|4*^|) . (p?5)) and ([26l in (|5^|) and rearranging terms, we arrive at the result 

/ 4 \ 3s 2 

^V 1 --^)'^ (53) 

on account of the fact that Q = s on the hypersurface E. This equation gives the relation between the heat flux, 
radial pressure, shear viscosity and the charge per unit length of the cylinder, over the hypersurface E. The result 
shows that for an uncharged radiating fluid without any shear viscosity, the radial pressure equals the heat flux all 
over the boundary of the collapsing cylinder. Equations (f5Tj) and ([53")) are generalizations over the results obtained 
earlier in (38| and (39|. 

The total luminosity of the collapsing matter visible to an observer at rest at infinity is 



fdM\ ( dM dt fdu^ ' ' 



V du 1 v \ dt dr \dr 



(54) 
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Differentiating (g9) with respect to t and using ([33]). (148]) . ((25)) and ([26]), we obtain 




= 4ttC" 



which, in view of (|53[) leads us to the expression 



3s 2 
87rC 4 



B 



(q + e) 





(55) 



(56) 



Thus the total luminosity of the collapsing matter as visible to a distant observer, depends on the energy flux associated 
with the collapse. Since Chiba 31| has demonstrated the absence of horizons in the cylindrical system, we expect the 
end state in this case to be a naked singularity. Therefore the collapse will terminate at the time when 




A. Dynamical Equations for the Collapse 



The dynamical equations for non-adiabatic charged anisotropic fluid with shear viscosity undergoing cylindrical 
collapse can be obtained from the Bianchi identities (T Q/3 + T^ m ^), p = for energy-momentum conservation. 
Using J2]), ©, and ([15]), we have 

(T- + K-^)^. - 4( A + i )-^( (l + P r + 2e --l w )-^( M + P ±+e+ ^ 

and 



1 f '2(q + e) ( B <j\ ss 1 

+ a^ + ^ + ^{b + c)-^bc-^ - < M > 

To discuss the dynamics of the collapsing system, it is customary to introduce the proper time derivative 

D T = (59, 
and the proper radial derivative Dr constructed from the circumference radius of a cylinder inside £ 

d* = 55f . w 

where 

R = C. (61) 
The fluid velocity for the corresponding collapse is given by 

U = D T {R) = (62) 

which must be negative to ensure collapse to occur. Defining new variable e = (note that e and e arc different 
quantities) and using ([H]) . we have 
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Consequently, Eq. (|51j) can be re-written as follows: 

4tt((7 + e) = 



(64) 



The time rate of variation of the total energy inside the collapsing cylinder is given by 



D T E' = -4wR 2 



Pr 



7! w 



1 



32irR 2 



U + e(q + e) 



3s 2 U 
2R 2 ' 



(65) 



In the case of collapse, since (U < 0), the coefficient of U inside the square brackets, will increase the C-energy of the 



cylinder if P r + e 



32nR 2 ' 



i.e. the effective radial pressure is greater than a certain value. The work done 



by the effective radial pressure leads to the increase of C-energy. The second term in the square brackets, due to the 
overall negative sign, describes the outflow of energy in the form of heat flux and radiation during the collapse. Since 
the collapsing cylinder contains the same species of the charges, the last term will decrease the energy of the system 

Q 2 

as plays the role of Coulomb force of repulsion and U < 0. 

The variation of energy between the adjacent coaxial cylinders inside the fluid is given by the expression 



D R E' = 4ttR 2 



U 



(<7 + e) 



I 



R 



Drs 



3s 2 
2R 2 ' 



(66) 



The first term on the right hand side gives the contribution of the energy density of the element of fluid inside a 
cylindrical shell, along with heat flux and radiation. Since U < 0, the factor ^{q + e) decreases the energy of the 
system during the collapse of the cylinder. In the remaining terms, the constant 1/8 comes from the definition of 
C-energy and the other term is the electromagnetic contribution. The C-energy of the cylinder at a given instant 
of time is then obtained by integrating (|66j) from the axis to the periphery of the cylinder. The acceleration of the 
collapsing matter inside the hypersurface £ is obtained using (TJJ}, (|2l>)) . (1521) and 



E' ) — 4vri? 
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Pr 



V3 



eA! 
AB 



5s 2 
2i? 3 ' 



(67) 



Substituting for 4- from Eg. (1571) into Eg. (155)) . we obtain the equivalent of Newton's second law of motion for the 
collapsing matter in the form 



4 \ / 4 

H + P r + 2e- -^WJ P>tU = -\^i + P r + 2e- -j=r)o 



4 \ 5s 2 



. 4{q + e)U n . . 1 / B C 



D R [P r + e - -^=w) + 2 (Pr - P± + e - 2V3 V o) 4 - ^jD r 



V3 



R nR 4 ' 



(68) 



which can be simplified as follows: 



(/x + P ref{ + 2e) D T U = - (n + P reff + 2e) 



D T (q + e) 



R 







(i- 


8)1 



DR(P rcff +e) + 2(P r ^ ff -P ±cff +e) 



1 



R nR 4 



Drs 



Here we assume that in general -4- (J| — ^) ^ 0. The left hand side of (|6"8")) represents force. The factor (/i + P r + 
2e — -^V 7 ) represents the inertial mass density, which gives the effect of dissipation but there is no contribution of 
the electric charge, nor of heat flux. The remaining term on the left hand side is acceleration. Thus, we can say that 
the dynamical system will evolve radially outward or inward according as DtU < or DtU > 0. The terms with a 
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negative contribution in (|68|) . favors the collapse while the other contribution prevents the collapse. If both of these 
terms cancel each other, then a condition of hydrostatic equilibrium will be encountered. 

The first term on the right hand side represents the gravitational force. The factor within the first square brackets 
shows the effects of specific length, effective radial pressure and the electric charge on the term (fx + P r + 2e— -^V 7 ) 
representing the active gravitational mass. The second term represents the contribution due to radiation and heat flux, 
which will leave the system (if there is an overall negative sign) through the outward radially directed streamlines. Thus 
it is in the same direction of pressure and would prevent the collapse. The third term has three main contributions: 
the first is the effective pressure gradient which is always negative, thereby preventing the collapse. The second is 
the local anisotropy of the fluid which will be negative for P r ^ ff < P. j .. fJ . in which case it will decrease the rate 
of collapse. The third is the electromagnetic field term. According to [2 if . the third term contributes negatively if 

> Drs. Under these conditions, the term in the third square brackets, with negative sign, contributes positively 
by reducing the attractive nature of the force appearing on the left hand side of this equation and hence this term 
will prevent the gravitational collapse. 



B. Solution of the Field Equations 

In their work, Di Prisco et al. [37} have derived the solutions for the shearfree and isotropic case of cylindrical 
collapse. Sharif and Abbas [46| have found analytical solutions for charged perfect fluid cylindrical gravitational 
collapse. The general solution in presence of dissipation, using the method of separation of variables, should be of the 
following form: 

A(t,r) = A (r), 

B(t,r) = B„(r)/i(t), (69) 
C(t,r) = C„(r)/ 2 (t), 

where Ao(r), Bo(r) and Co(r) are solutions of a static fluid having /j,q as the energy density and p r o and p±o as 
the radial and tangential pressure. These solutions for such a general case, are highly complicated and hence we 
concentrate on a simplified choice for the facility of calculation. For that we extend the solutions assumed in [l8| for 
the spherically symmetric distribution to the case of dissipative cylindrical collapse with the choice of f± (t) = 1 and 
h{t) = f(t), in Eq. (JBSJ) above, i.e. 

A(t,r) = A {r), 

B(t,r) = flo(r), (70) 
C{t,r) = C (r)f(t). 

On account of equation ((33 j) in the junction condition, the function C(t,rs) represents the luminosity radius of the 
collapsing matter as seen by a distant observer. The expression f) 12|) for the shear scalar becomes 



1 / 



V3Ao /' 

The field equations (|M|) to ([27)) are reduced to 



(71) 



87r(/i + e) = 87r/xo + ^-^rp (72) 



M? + e) = - r ^- ? f^-4a , J (78) 



A Bof \C A c 




8^=8^0-^ + — -^ (75) 



11 



where 



Mo 



2C ' 




(2B' 


C 


C 


V Bo 



9k 

Co 



(76) 




2A' C'o 



(77) 



- ■ 


A)C y 






C A 


-BoCo\ 




c y 


-B Co / 



and 

RttP, „ = , 

S Q \^4o Co -B 

Equations (l72j) to (|75|) represents the static anisotropic fluid configuration in the limit f(t) — > 1. Substituting (|71[) . 
(|73|) and (|74]l into ([53]) and assuming that P r o = 0, we obtain the following differential equation: 

2// + ./ 2 + a//-^ = 0, (79) 



where a and b are functions of r, i.e 

n = , 

-Bo V Co Ao 



2) (80) 



and 



A 



2 ,,2 



Co 4 



(81) 



Equation (f79f can be solved numerically assuming suitable boundary conditions. Although the first three terms on 
the left hand side are similar to Eq. (56) in [l8j], but this case is clearly distinct due to the difference in geometry and 
the nature of collapsing matter. Thus the parameters a and b are different in their content. The parameter b in our 
case, depends on the charge enclosed inside the cylinder and hence may be positive, negative or zero. 
The solution is sought under the following boundary conditions: 

1. Initially the system represents the static configuration at t — > — oo, when f{t) — > and f(t) — > 1. 

2. Further, if possible, as t — > 0, f(t) — > 0. This means that the luminosity radius C(t, r^) has the value Co(r^) at 
the beginning of the collapse and vanishes as the final fate is approached. 

Using the second boundary condition in (|79|) . we find that f(t) becomes infinite as t — > 0. Thus the luminosity will 
remain fairly constant during the collapse and will suddenly go to zero value at the end. Moreover, applying the 
initial condition we find that the second order derivative of f(t) i.e. f(t) depends on the nature of the charge enclosed 
inside the cylinder at t — > — oo. Therefore we can divide our subsequent analysis into two parts: 



When the fluid is uncharged: In this case the last term on the left hand side of ([79| vanishes irrespective of 
the instant of time. Thus ([79"]) is now reduced to the form 2/ / + f 2 + aff = 0. It is a second order, reducible 
differential equation of the Liouville type and can be solved using Maple or Mathematica program. The solution 
is a decreasing function of time, with the rate of decrease being dependent on the parameter a and is given by 

f(t) — (^y 1 — 322. e - Qt / 2 ) 2 / 3 j where C\ and C2 are arbitrary constants. Expanding binomially, neglecting the 
higher order terms and choosing suitable values for C\ and C2, we find that f(t) ~ 1 — exp(—at/2). 

When the fluid is charged: In the case of a charged fluid, the parameter b is non-zero, and the solution 
is not straightforward. If this charge is negative, then the second derivative will be negative at t — > —00 and 
hence f(t) will be a decreasing function of time as t — > 0. When the charge inside the cylinder is positive, 
f(t) will increase as t — ¥ 0. To achieve a numerical solution, we modify our trial solution in the form A(t,r) = 
Ao(r),B(t,r) = B (r) and C(t,r) — A (r)f(t). This corresponds to the case for which (q + e) = 0, although 
q and e are separately non-zero. This feature arises due to the special type of geometry which is involved in 
this case and is not a general feature. The evolution equation is now: 2ff + f 2 — i = 0, which can be solved 
numerically by assuming sample values for b, but for that we have to discard the second boundary condition, 
i.e., we now have at t — > 0, f(t) ^ 0. 



Sample plots are shown in the following figures. 



12 




FIG. 1: Diagram showing the plot of f(t) vs t for the uncharged case with increasing values of a in the left figure, the charged 
case with b < in the middle figure and with b > in the right figure. The magnitude of b is the same in the last two figures. 



V. SUMMARY 



We have formulated, by the help of General Relativity, the field equations for the case of dissipative cylindrical 
collapse in presence of heat flow, free-streaming radiation, and shear viscosity and have obtained a few results. We 
have derived the relation between the expansion 9, the shear a and the energy flowing out of the cylinder in the 
form of heat flux q and free-streaming radiation. By employing the Darmois-Israel junction condition for the smooth 
matching of interior and exterior spacetimes at the boundary S, we have verified the relation between the specific 
energy of the cylinder due to the electromagnetic field and the mass of the collapsing matter. The total luminosity of 
the collapsing matter as visible to a distant observer, depends on the energy flux associated with the collapse. This 
energy flux over the hypersurface S bounding the cylinder, is dependent on the effective radial pressure and the charge 
per unit length of the cylinder. The dynamical equations for the collapse is derived from the Bianchi identities with 
the help of Misner-Sharp formalism for this non-adiabatic, anisotropic and dissipative fluid and the equation for the 
effective Newton's second law of motion is constructed. Finally, we have presented the solution to the field equations 
for the given matter distribution undergoing cylindrical collapse and have interpreted its significance from a physical 
point of view. 
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